Mon. Not. R. Astron. Soc. 000, 000-000 (1994) Printed 15 September 1995 



Accretion discs with accreting coronae in AGN — 
II. Nuclear wind. 

Hans J. Witt/ Bozena Czerny^ and Piotr T. Zycki^ 

^ Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 144^^ Potsdam, Germany, e-mail:hwitt@aip.de 
^ Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland, e-mail:(bcz,ptz)&camk. edii.pl 



ABSTRACT 

We study an accretion disc with a hot continuous corona. We assume that the corona 
itself accretes and therefore it is powered directly by the release of the gravitational 
energy and cooled by radiative interaction with the disc. We consider the vertical 
structure of such a corona and show that the radial infall is accompanied by strong- 
vertical outflow. Such a model has two important consequences: (i) at a given radius 
the corona forms only for accretion rate larger than the limiting value and the fraction 
of energy dissipated in the corona decreases with increasing accretion rate, and (ii) 
the disc spectra are significantly softer in the optical/UV band in comparison with 
the predictions of standard accretion discs due to the mass loss and the decrease of 
internal dissipation in the disc. Both trends correspond well to the mean spectra of 
radio quiet AGN and observed luminosity states in galactic black hole candidates. 

Key words: accretion, accretion discs - galaxies: active - galaxies: Seyfert - X-rays: 
galaxies 



1 INTRODUCTION 

Hard X-ray emission in Active Galactic Nuclei results most 
probably from the Compton upscatter of soft photons by a 
hot optically thin plasma but the place of origin of this emis- 
sion is still under debate despite the amount of observational 
data accumulated (e.g. Bregman 1994, Czerny 1994). 

A number of models were suggested. The model of 
a hot optically thin continuous corona above an optically 
thick accretion disc first discussed by Liang & Price (1977) 
is particularly attractive and significant progress has been 
recently made along this line starting mostly from the pa- 
per of Haardt & Maraschi (1991), followed by Nakamura 
& Osaki (1993), Kusunose & Mineshige (1994), Svensson 
& Zdziarski 1994, Zycki, CoUin-Soun'rin & Czerny (1995; 
hereafter Paper I). It was shown that the dissipation of a 
significant fraction of the gravitational energy in the corona 
tends to stabilize the disc (lonson & Kuperus 1984, Svens- 
son & Zdziarski 1994). A number of papers were devoted to 
the stabilizing role of advcction in the disc surrounded by a 
corona (e.g. Chen et al. 1995). The advection-dominated 
models were successfully applied to some galactic X-ray 
sources (Abramowicz et al. 1995). 

Most of the previous papers were based on the assump- 
tion that the fraction of the energy generated in the corona 
is a free parameter of the model and does not depend on 
the accretion rate. This value was assumed to be constant 
also during time evolution computations (Abramowicz et al. 
1995). 



However, conclusions from the model crucially depend 
on this assumption (Paper I). A coupling between the frac- 
tion of the energy dissipated in the corona and the accretion 
rate may lead to different picture. 

Observations suggest that such a coupling docs exist. 
Direct evidence is available only for some galactic X-ray 
sources, namely for X-ray novae. They show rapid change 
of the spectrum at the end of an outburst from soft X-ray 
dominated spectrum to the domination by hard X-ray power 
law (e.g. Miyamoto et al. 1993). This suggests that when 
the accretion rate drops to ~ 1% of the initial value, a tran- 
sition occurs from a disc-dominated solution to a corona- 
dominated solution (Miyamoto et al. 1993, van der Klis 
1994). The initial value is probably close to the Edding- 
ton limit. 

There are no such clear indications for a similar trend 
in Seyfert galaxies or quasars; however, the hardening of the 
UV spectrum of Seyfert galaxies with the increase of the lu- 
minosity may indirectly support similar correlation between 
the fraction of energy dissipated in the corona and the ac- 
cretion rate. 

Therefore a method to predict theoretically the fraction 
of energy liberated in the corona is of extreme importance 
as its predictions can be tested against observations thus 
allowing for deeper insight into the very phenomenon. 

Such a method was outlined in the Paper I. However, 
that paper was based on the assumption that the corona 
is in a hydrostatic equilibrium which was not justified for 
higher accretion rates so conclusions from the model could 
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not have been drawn. 

In this paper we study a corona above an accretion disc 
allowing both for the vertical supersonic outflow as well as 
the radial advection term, with the description of the cou- 
pling between the disc and the corona as in Paper I. The 
description of the vertical structure of the corona including 
transonic vertical outflow is given in Section 2. In Section 
3 we show the properties of the corona including the de- 
pendence of the calculated fraction of the energy liberated 
in the corona on accretion rate and viscosity, as well as we 
estimate the mass loss from the disc and the radial depen- 
dence of the disc accretion rate. We also compute exemplary 
disc/corona spectra in the optical/UV/X-ray range taking 
into account the mass loss from the disc. In Section 4 we 
discuss the importance of the radial advection. Conclusions 
are given in Section 5. 

2 DYNAMICS OF THE CORONAL OUTFLOW 
2.1 Equations of vertical motion 

Following Begelman, McKee & Shields (1983) and Meyer & 
Meyer-Hofmeister (1994) we derive a one-dimensional model 
for the radial accretion and outflow process in a cylinder 
along the «-axis. We assume that the gas is isothermal and 
the outflow is driven by radiation pressure from the disc. 

Starting with the continuity equation of mass flow in 
cylindrical coordinates we have 
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where Vr and Vz are the radial and vertical components of 
the velocity. The radial flow velocity can be approximated 
by 
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Here Cs denotes the isothermal sound velocity (cg = P/p), 
Qk is the Keplerian angular velocity, G is the gravitational 
constant, M the mass of the black hole and a is the viscosity 
parameter. 

We replace the sidewise loss term of the vertical flow by 
(cf. Meyer & Meyer-Hofmeister 1994) 
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The vertical flow is not restricted to the z-axis and the 
flow lines diverge while the geometry changes asymptotically 
from the cylindrical to spherical for z > r. This divergence 
is essential for the transonic outflow to develop and cannot 
therefore be neglected. To account for this effect, we replace 
the vertical mass flux with the mass flux in the diverging 
tube (e.g. Begelman et al. 1983) 
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The introduced assumptions allow us to rewrite the conti- 
nuity equation (1) in the following form 

d{pVz) _ 2pVrr _ 2zpVz , s 

dz ~ r'' + z'' r'' + z''' ^ ' 

where m{z) = pVz now describes the vertical accretion rate 
per unit area. It depends only on z. 



Second differential equation describes the motion of the 
gas in the corona in the vertical direction. 
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where P is the gas pressure, the last term is due to the grav- 
ity of the central object (we neglect any self-gravity effects 
in the disc) and the radiation pressure gradient in the corona 
is given by 

dPiad _ FsottKesP 

dz ~ c ' ^> 
where Fgoft is the total flux emitted by the disc resulting 
from internal viscous generation and tliermalisation of the 
external irradiation, and Kes = 0.34cm^/g is the electron 
scattering opacity. 

Replacing now all terms pVz by niz and using the rela- 
tion 

dvz diriz rhz dP 



dz 



dz 



P dz ' 



we obtain two differential equation for mz{z) and P{z) 
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where T is the corona temperature (equal to the ion temper- 
ature T in the case of two-tempcraturc plasma; sec Section 
2.3). These equations describe the equilibrium between ac- 
cretion and outflow due to radiation pressure coming from 
inverse Compton scattering. 

These two differential equations can be partly solved 
analytically as shown in the Appendix A. 

The values of the temperature of the corona and the soft 
(disc) flux has to be determined from the energy balance (see 
Section 2.3). 

2.2 Basis of the corona and the sonic point 

The boundary conditions for the differential equations (9) 
and (10) are given by the basis of the flow pattern being the 
surface of an accretion disc and by the condition that the 
outflow has to be transonic. 

The height of the disc or the basis of the corona is de- 
fined by the equilibrium of the radiation pressure and the 
gravitational force in vertical direction. 
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The equation has two solutions for z if Fsoft is smaller than 
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For the limiting value the two solutions merge to one, z = 
r/V2. If the flux Fsoit exceeds this limit a disc may not 
form because the system is not in hydrostatic equilibrium. 
This limit corresponds to the Eddington luminosity but the 
numerical coefficient is slightly different due to the flat (and 
not spherical) geometry. 
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The physically relevant solution for the height of the 
disc, zo, is taken from the range 0<zo< r/V2. The other 
solution is dynamically unstable. 

The condition of the regular behaviour of the outflow 
at infinity is replaced, as usual, by the condition that the 
proper solution of the two differential equations (9) and (fO) 
has to go through a sonic point (cf. also Bondi 1952). A sonic 
point is a removable singularity where both the nominator 
and the denominator vanish simultaneously. We therefore 
obtain the condition Pgonic = "^sonic^^/^^H- Inserting this 
condition in the nominator of equation (10) yields 
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This means that wc can determine the position of the sonic 
point without knowing the explicit solution of P{z), which 
must go through the sonic point. The sonic point Zsonic 
depends only on some external quantities like the corona 
temperature T, viscosity parameter a, radial coordinate r, 
etc. 

Mathematically speaking equation (13) represents only 
an extension of equation (11) in which two positive terms are 
added (consider that < by definition). Equation (13) 
has not more than two solutions, like equation (11). Since 
the presence of the two terms on the left side is equivalent to 
an increase of -Fgoft in equation (11), it follows that Zsonic > 
Zo- In addition we can notice that the sonic point must 
disappear if the temperature becomes too large. 



2.3 Energy budget 

Since we assume that the corona is isothermal (i.e. there 
is no temperature gradient in the vertical direction) our de- 
scription of the energy generation and cooling basically (but 
not fully) follows Paper I. 

We assume that both the corona and the disc dissipate 
the potential energy of the accreting gas. We therefore adopt 
the a-viscosity description for both layers. We assume a two- 
temperature plasma of ions and electrons. In a steady state 
the heating and cooling rates must be equal. The dissipation 
energy of the ions is transported via Coulomb interaction to 
the electrons. The heating of ions arises from the viscous dis- 
sipation of the gravitational energy, while the cooling is con- 
ducted by Coulomb electron-ion exchange. Electrons heated 
by the last process cool by Compton scattering off soft pho- 
tons from the disc. The soft photon flux comes partly from 
the energy generated in the disc and partly from the fraction 
of hard X-ray photons downscattered by the corona, as in 
the basic paper of Haardt & Maraschi (1991). 

The basis of the corona is determined by the condition 
that the ionization parameter, H, (Krolik, McKee & Tarter 
1981) marginally allows the medium to reach the inverse 
Compton temperature (i.e. the atomic processes are negligi- 
ble). 

These conditions can be formulated in such a way that 
the structure and properties of the cool disc need not to be 
known explicitly. 



The total generated flux of the model is given by the 
global parameters the black hole mass M, the global accre- 
tion rate M, the viscosity parameter of the corona a and the 
radial coordinate r where the one-zone model of the disc is 
considered. 

The total flux is then given by 
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where /(r) represents the boundary condition at the margi- 
nally stable orbit 

/(r) = l-(3i?schw/r)''" (14a) 

in the Newtonian approximation. 

For the flux generated in the disc we simply set 

FA = ^Ftct=aFd + Fc), (15) 

where the fraction of energy dissipated in the disc, ^, has to 
be calculated by iteration using the total flux and the flux 

generated in the corona by a-viscosity. The flux generated 
in the corona by dissipation is given by 
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The cooling of ions in the corona by the electron-ion 
Coulomb interaction is described by the following equation 
(Shapiro, Lightman &; Eardley 1976) 
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where Ti, Te are the ion and electron temperatures, and 



2.44 X lO^V^T^'^ 



In A; with In A « 20 



is the electron-ion coupling rate. 

For the soft flux fsoft we have the relation 

fsoft(«o) = + »?-Fc(l - a). 



(17a) 



(18) 



where r] is the fraction of the coronal flux directed towards 
the disc and a the X-ray albedo. We assume for our model 
r] = 0.5 and a = 0. The soft flux results in Compton cooling 
of hot electrons which is balanced by heat generation. 
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where y is the Compton parameter for optically thin 
medium, 

y = Tes + J : (20). 

and Tes denotes the scattering optical thickness of the 
corona, 



oo 



dz. 



(21) 



Finally we require that the value of the ionization pa- 
rameter H at the basis of the corona corresponds to marginal 
thermal stability: 

" = "°llO^J =2^- ^^^^ 
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Here Pq = P{zo) denotes the pressure at the basis of the 
corona and Ho = 0.65 is the ionization parameter for in- 
verse Compton temperature 10^ K. We take into account 
that the critical value of that parameter depends on the in- 
verse Compton temperature (Begelman et al. 1983). 

Therefore in our present approach we made two impor- 
tant changes with respect to the thermal conditions formu- 
lated in Paper I: 

(i) we define H using corona gas pressure instead of the disc 
gas pressure 

(ii) we scale H with the electron temperature in the corona. 
Both changes lead to essential changes of the character 

of the model. 

Equations (9) and (10) with proper boundary condi- 
tions, supplemented by equations (14)-(22), allow determi- 
nation of the structure and the dynamics of the hot phase 
including the fraction of energy liberated in the corona. The 
solution at a given radius is parameterized by the mass of 
the central object, the accretion rate and the viscosity pa- 
rameter Q. 



2.4 Cold disc structure and the stability of 
corona/disc solution 

Although the solution is fully determined by the equations 
for the corona we may be also interested in calculating 
the properties of the cold layer in equilibrium with such a 
corona. The structure of the disc can be uniquely calculated 
by integrating equations of the disc vertical structure with 
the surface boundary conditions imposed by the corona so- 
lution (Czerny, Zycki & CoUin-Souffrin 1995). Simpler esti- 
mates based on the vertically averaged disc structure can be 
obtained but they are less reliable as the departure of opac- 
ities from pure electron scattering and the pressure profile 
influence the results. 

Computations of the disc structure allow to calculate 
the total (i.e. disc plus corona) surface density as a function 
of the accretion rate for our model. This function is closely 
related to thermal stability of disc/corona solutions (Bath 
& Pringle 1982). 

A coronal solution (i.e. the values of ^ and Po) gives the 
boundary conditions necessary to solve the differential equa- 
tions of the disc vertical structure. We do the computations 
(see Section 3.5) assuming the aPtot description of viscos- 
ity. In principle, the viscosity coefficient in the disc may be 
different from the viscosity coefficient in the corona as the 
thermal conditions are actually very different in these two 
parts of accretion flow. However, as two values of a would 
only lead to quantitative, and not qualitative, differences 
from a single value model we assume the same value for the 
two coefficients in the present computations. 



3 RESULTS FOR THE CORONA WITHOUT 
RADIAL ADVECTION 

3.1 Topology of solutions; transonic curves 

We first discuss the topology of solutions of the two differen- 
tial equations (9) - (10) and the conditions of the existence 
of transonic curves without referring to algebraic or integral 
equations describing the thermal balance. 



Figure 1. Two examples of topology of transonic solutions. The 
parameters are indicated in each figure. In the top panel (a) the 
correct transonic solution marked by a thick curve determines 
uniquely the initial value of Po at the basis of the corona. In the 
bottom panel (b) a case is shown where the transonic solution is 
not connected to the basis of the corona. 

Detailed study of the equation (13) determining the po- 
sition of the sonic point shows that it may have two solu- 
tions. The two solutions merge with increasing coronal tem- 
perature (the ion temperature, Ti, in our two-temperature 
approximation), or the soft flux from the disc Fsoft, and then 
the solution disappears. 

For T ^ the solutions of equation (13) converge 
towards the two solutions of equation (11), zo and zi « 
(r/zo — 3/2)"'/^. The solutions of equation (13) are confined 
between zo and zi. Therefore, if the soft flux is increased, 
the sonic point disappears before the system reaches the Ed- 
dington limit as defined for spherical symmetry. The merg- 
ing point yields an upper limit of the temperature in the 
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Figure 2. The positions of ttie two sonic points (see equation 
13) as functions of the corona temperature T; for three values of 
5 (or Fgoft; see equations 14, 15 and 18). 

corona which cannot be exceeded, 
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The formula is derived for a = and small zq. The maxi- 
mum possible temperature in the corona is smaller for higher 
viscosity parameter a. 

The solution close to the disc surface is positioned at 
zo < z < Emerge, where the merging point cannot exceed 



the limit 



< 0.835r. The corresponding singularity 



is of the focus type. There is no transonic solution pass- 
ing through such a singularity. Physically it is related to 
the fact that at this position the gravitational potential still 
increases outwards so there is no acceleration of the out- 



flow. The other solution is located at Zn 



< z < zi where 



the gravitational potential decreases outwards so the veloc- 
ity becomes transonic. This type of singularity is called a 
saddle and it allows to cross from subsonic to supersonic 
motion. The coexistence of the two singular points is a typ- 
ical property of a flow with a potential barrier, e.g. in the 
case of disc accretion in the vicinity of the marginally stable 
orbit (Muchotrzeb 1983). 

Two examples of the topology of the solutions of equa- 
tions (9)-(10) are shown in Fig. 1. These figures were de- 
rived by varying the integrating constant (sec for details 
Appendix B), i.e. a number of different initial values of the 
pressure was chosen. In the top panel (a) we show an ex- 
ample of the existence of the disc/corona solution. There 
is a curve which starts at the disc surface at high pressure 
and low velocity values, and passes through the sonic point 
at 2; « 17.7 achieving, with increasing z, low pressure and 
high velocity values. In the bottom panel (b) we show a sit- 
uation where the solution which is going through the sonic 
point is not connected to the basis of the corona. For such 



Figure 3. Fraction of the energy liberated in the disc as a func- 
tion of accretion rate; the corona dissipates a fraction 1 — 5. Lower 
limit for accretion rate corresponds to g = whilst the upper limit 
for accretion rate and 5 is due to the disappearance of transonic 
solutions. The curves are labeled by the value of the viscosity pa- 
rameter a. The dashed curves show the analytical approximation 
given by equation (25) for a — 0.1 and 1. Other parameters are: 
M = IO^Mq, r = lOiischw 

parameter space there is no uniquely physical solution. For 
the following discussion it is important to understand that 
the family of physical solutions is limited by the fact that 
either the two sonic points converge or the transonic curve 
cannot be accessed from the disc surface. 

The convergence of the two sonic points is illustrated in 
Fig. 2. 

3.2 Fraction of energy liberated in the corona 

When the results from the previous section are supple- 
mented by the equations of energy budget (see Section 3.3), 
the soft flux or, equivalently, £, are no longer free parameters. 

In Figure 3 we show the dependence of ^ on the accre- 
tion rate. The computations were done for central black hole 
mass 10* Mq, at a radius of 10 -Rschw The curves correspond 
to different values of the viscosity parameter a ranging from 
0.1 to 1. 

The most striking differences between the relation ^ vs. 
m obtained here and in Paper I are: (i) the existence of 
one solution here instead of two in Paper I and (ii) the ba- 
sic solution of Paper I predicted the increase of the relative 
importance of the corona with accretion rate whilst our so- 
lution predicts the opposite trend, (iii) the solutions exist for 
much broader, albeit still limited, range of accretion rates. 

These differences result mostly from the two changes 
introduced to the consideration of the energy budget and not 
from the inclusion of the dynamical terms into the corona 
equations. 

First, we replaced the gas pressure at the disc surface 
with the corona pressure at the basis in relation (22). Disc- 
dominated models obtained by Paper I had ^ ~ 1 and the 
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Figure 4. The dependencies of corona parameters: the ion temperature, Ti, the electron temperature, Te, the electron scattering optical 
depth, Tes, and the Compton parameter, y, on the accretion rate, M, for a number of values of the viscosity parameter a. Model 
parameters: M — 10* Mq, r — lOiigchw 



disc surface pressure was strongly dominated by the radia- 
tion pressure so the total pressure equilibrium did not result 
in gas pressure equilibrium at the disc/corona boundary. 
Our new definition of the disc/boundary layer seem to be 
more correct as we use tlie location of tlie upper bend on tlie 
T vs. H curve to define the transition (Krolilc et al. 1981; see 
also Paper I). However, the problem indicates the necessity 
of the more detailed computations of the transition zone. 

Using our present definition of H witli some other sim- 
plifications (Ti ^ Te, He <S r, 4kTe/meC^ <C 1), we obtain 
an approximate relation for H = Ho, 

M = A{l-0-' (I'lf^) ' (24) 
where 

A = 2.5xlO''El^/'a-'Q^' fir)-' [g/s]. 

This equation shows that the fraction of the energy gener- 
ated in the corona increases with the accretion rate. How- 
ever, the assumption H = const leads to relatively narrow 
range in m with interesting solutions. These solutions exist 
only for relatively low values of 5o (~ 0.01); otherwise the 
accretion rate exceeds the Eddington rate. 

Much broader range of discs with relatively strong coro- 
nae exist if the dependence of E on the coronal temperature 
is included, as we do in the present model (see equation 22). 
An approximate solution can be derived also in this case: 

M = A{l-0-'(ln^^ , (25) 

where 

^ = 8.3 X 10'' a-' Q^' fir)-' [g/s] 



Although these solutions approximate relatively well the ex- 
act solutions for intermediate values of ^, they misrepresent 
them for both very high and very low accretion rates. For 
very low accretions tlio correction terms due to tlic liigli 
value of the electron temperature are important, since the 
electron temperature varies as Te ~ m~'^'^ (see Appendix 
C). 

In the case of high accretion rates the dynamical terms 
become important. However, tlieir inclusion in the corona 
structure does not extend tlie existence of solutions into 
much higher accretion rates than that allowed by tlie liydro- 
static equilibrium models (see Paper I). The physical mech- 
anism limiting the accretion rate is however different from 
the one discussed by Paper I in the case of liydrostatic equi- 
librium solutions. Now the two sonic points converge (see 
Section 3.1) with increasing m because both T and Fgoft 
increase with m. There are no stationary transonic solu- 
tions above the merging point i.e. for higher m, and this 
conclusion is not ciianged if tlie advection term is taken into 
account (see Section 4) The limit depends on the viscosity 
parameter a, as can be seen from Fig. 3. 

With decreasing accretion rate the relative importance 
of the corona increases. The limit ^ = 0, i.e. all energy 
generated in the corona (in this case the disc radiates only 
at the expense of irradiation by corona), is reached for finite 
value of the accretion rate. For radial distance 10 i?schw and 
M = 10** M0 this limit is ~ 8 x 10^^ g/s. For lower accretion 
rates there are no solutions for a disc with a corona. 

As tlie fraction of energy liberated in tlie corona de- 
pends explicitly also on the current disc radius, this fraction 
is a function of a radial distance for a given accretion rate. 
We sliow tlie corresponding plot (Fig. 6) under tlie assump- 
tion that the accretion rate is constant, i.e. the vertical losses 
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are relatively unimportant. The corona covers only a lim- 
ited fraction of a disc (r < rmax) which increases with the 
accretion rate. It develops mostly near rmax and its relative 
importance decreases inwards. 

3.3 Properties of the corona 

In Fig. 4 we illustrate the properties of the corona. The 
computations are made at a radius 10 i?schw, assuming M = 
10** Mq. 

The optical thickness of the corona is moderate (~ 
0.1 — 0.2) so the assumption of optically thin medium is 
well satisfied. The electron temperature is never very high 
so there is no significant pair creation process. The ion tem- 
perature is mostly higher than the electron temperature al- 
though the difference between them is lower for low accretion 
rate. As a result, the value of the Compton parameter y is 
never very high; it is close to 1 only for accretion rate close 
to the minimum value and is decreases from 0.6 to 0.05 for 
accretion rates increasing from 10^^ to almost lO^^g/s. 

Approximate analytical solutions illustrating the depen- 
dence of the coronal properties on the disc radius are given 
in Appendix C. Those formulae are not very accurate be- 
cause at higher accretion rate the departure from hydro- 
static equilibrium becomes essential and at low accretion 
rate the nonlinear terms in equations(17) and (20) are im- 
portant. However, as the dynamical effects are not impor- 
tant at low accretion rates approximate numerical results 
can be found without the necessity to solve numerically the 
transonic outflow problem (see Section 3.4 and Appendix 
D). 

3.4 Vertical vs. radial flow and global solutions 

The relative importance of the vertical outflow can be 
conveniently measured by the asymptotic value, mj°, of the 
expression mz{l + z^ /r^) normalized to a logarithmic deriva- 
tive of the mass flux at a given radius, i.e. 



dlnM 47rr^m° 



dlni 



M 



(26) 



Numerically computed values at 10-Rschw as a function of 
accretion rate are shown in Figure 5. We see that the mass 
loss strongly depends on the accretion rate - it increases sig- 
nificantly with the luminosity assuming constant viscosity 
parameter a. What is important, it also strongly depends 
on the value of the viscosity parameter a. For a close to 1 
mass loss is never of any importance whilst for a of order of 
0.1 it may cause the outflow of the major (dominating) frac- 
tion of the mass flux. As the value of viscosity is not known 
the detailed predictions of the mass loss are impossible along 
this line. We can only show a few examples of the possible 
consequences of the outflow and, if the overall detailed pre- 
dictions of the model are determined and confronted with 
the data for a range of values of a, the viscosity parameter 
can be constrained independent from other estimates (e.g. 
Siemiginowska & Czerny 1989). 

We now consider in some detail the case of a = 0.1. The 
mass loss from the corona is clearly important and the global 
solution for the disc cannot be obtained assuming M = const 
in the region where a corona develops. The correct solution 
has to be obtained numerically by computing the mass loss 
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Figure 5. The dependence of the vertical outflow ratio as a 
function of accretion rate (lower axis) or the fraction of energy 
liberated in the corona (1 — 5) (upper axis). Curves are labelled by 
the viscosity parameter a. Other parameters are : M = 10* Mq, 
r = lOiischw 

at every radius and then the disc/corona structure for the 
reduced accretion rate. 

As the construction of the exact solutions at every ra- 
dius and a range of accretion rates (as now the local accre- 
tion rate is not known a priori) is time-consuming we use for 
that purpose half-numerical and half-analytical approach to 
estimate the local mass loss rate. 

We compute first the corona parameters neglecting dy- 
namical terms (but treating carefully the equation of hydro- 
static equilibrium); we also include all nonlinear terms (see 
equations 17 and 20) as the mass loss in very sensitive to the 
obtained ion temperature. The set of equations (D4)-(D7) 
is given in the Appendix D. These equations yields us the 
approximate global solutions ^,Ti,Te and Po as functions of 
the parameters M, M, a and r. 

Using now equations (11) and (13) we can find the basis 
of the corona zo and the sonic point -Zsonic. Finally with 
equation (A6) we can estimate the pressure at the sonic 
point Psonic which subsequently yields msonic. 

Now we can use very good approximation. 



dlnM 
dlnr 
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which gives us the mass flux at a given radius r. 

This local value of the mass loss as well as the deter- 
mination of the local fraction of the energy liberated in the 
disc is incorporated into the code calculating the spectrum 
of a Keplerian disc radiating locally as a black body. Com- 
putations are done inwards which is equivalent to numerical 
integration of equation (26) starting from a given external 
accretion rate. 

As the fraction of energy generated in the corona in- 
creases for decreasing accretion rate at a given radius but it 
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Figure 6. Fraction of the energy liberated in the disc, 5, as 
a function of radius, r, for two values of the external accretion 
rate, 0.001 and 0.01 in units of 10^^ g/s. The dashed curves repre- 
sent the solution with vertical mass loss neglected whilst the solid 
curves show the solution with the effect of inward-decreasing ac- 
cretion rate taken into account (the solid and dashed curves over- 
lap for the lower value of M). The viscosity parameter a is equal 
0.1. 

also increases for a constant accretion rate and decreasing 
radius (see equation 25) we have two competing trends. The 
net effect is presented in Fig. 6. We see that the relative im- 
portance of corona decreases inwards although more slowly 
than in the case of negligible mass loss from the disc surface. 

The presented results, however, should be treated only 
as an indication of a trend, and not as strictly quantitative 
predictions. Further work towards better approximation of 
this two-dimensional flow is necessary. 

3.5 Thermal stability of the disc/corona system 

The total surface density of the disc/corona system resulting 
from integration of the vertical structure of the cool part of 
the flow (see Section 2.4) is shown in Fig. 7. 

The slope of the presented curve is negative for higher 
accretion rates. Therefore, although the presence of the 
corona suppresses the well known thermal instability of ra- 
diation pressure dominated discs, due to the irradiation (e.g. 
Czerny, Czerny &; Grindlay 1986) and the decrease of energy 
liberation (e.g. Svensson & Zdziarski 1994), it does not re- 
move the instability for higher accretion rates in our models. 

The problem is not a particular property of our model, 
but a general difficulty to accommodate two opposite re- 
quirements into the model. As the data clearly suggest (see 
Section 3.6.1) the UV bump should be more profound at 
higher accretion rates, i.e. the corona should be weaker at 
higher accretion rates. But weaker corona cannot stabilize 
the disc. Therefore, at least within the frame of the aPtot 
description of the disc interior, either we can, at least in 
principle, construct models which produce stable discs with 
stronger corona for higher accretion rates or unstable discs 



Figure 7. Relation: total surface density of the coupled 
disc/corona system, E, vs. accretion rate, M. Positive slope of the 
curve corresponds to thermally stable solutions whilst a negative 
slope indicates thermal instability at a given stationary accretion 
rate. The viscosity parameter a — 0.1. Other parameters are : 
M = 108MQ,r = 10iischw 

with weaker corona for higher accretion rates. The first 
type of models is consistent with observed spectral trends 
but opens the question of the development of thermal insta- 
bilities within the disc, particularly the problem of the final 
state. The second type of models give stable models but with 
spectral trends inconsistent with the data. This kind of dif- 
ficulties led to suggestions of the viscosity parameterization 
based on the 

Q^-f" gas 

model wdth variable (basically bimodal) 
Q (Chen & Taam 1994) leading to two stable branches and 
an unstable branch at intermediate accretion rates. 

3.6 Relevance to observations 

3.6.1 Big hump versus hard X-ray luminosity 

The observed multi-wavelength variability in many Seyfert 
galaxies and radio quiet quasars is characterized by harden- 
ing of the UV spectrum with an increase of brightness (e.g. 
Ulrich 1989) and mostly constant hard X-ray slope, partic- 
ularly if corrections for the warm absorber and/or reflection 
from cold gas is taken into account (see Mushotzky, Done 
&; Pounds 1993). Extensive monitoring campaigns for a few 
sources suggest that two distinctive kinds of variability are 
observed. Short term variability (days - months) is gen- 
erally characterized by UV flux being well correlated with 
the 2-10 keV flux (see Clavel et al. 1992 for NGC 5548. 
However, sudden increase of the UV or soft X-ray flux is 
sometimes observed without the corresponding increase in 
hard X-rays (for NGC 5548 see Clavel et al. 1992; Done et 
al. 1995). 

The first kind of variability is most probably due to 
some instabilities in the corona dominated fiow assuming 

that the UV radiation originates mostly from reprocessing 
of X-rays (e.g. Rokaki, CoUin-Souffrin & Magnan 1993). The 



second kind of behaviour reflects rather sudden increase of 
accretion rate in the innermost part of accretion disc; the 
timescale of 1 week for NGC 5548 is consistent with viscous 
timescale of the innermost parts of the disc for viscosity 
parameter a ~ 0.1, assuming the disc thickness to the ra- 
dial distance ratio about 0.005 (Siemiginowska, Czerny & 
Kostyunin 1995), the value of the central mass 3.7 X IO^'Mq 
(Clavel et al. 1992) and using the thermal time scale as de- 
rived by Siemiginowska &; Czerny (1989). 

Our model correctly reproduces the trends accompany- 
ing the increase of accretion rate. Particularly if the accre- 
tion rate is close to the minimum value, then minor changes 
in accretion rate can drive major spectral changes. 

3.6.2 UV spectrum 

The existence of the corona modifies the shape of the op- 
tically thick emission through the decrease of the energy 
dissipation in the disc as well as through Compton scatter- 
ing of the outgoing photons. In this section we concentrate 
on the first effect. 

According to our model, the optically thick emission 
(i.e. soft flux; see Section 2.3) is determined by the local 
value of ^ and M. These quantities are determined by equa- 
tions (27) and (CI) and the mass loss from the disc is calcu- 
lated as described in Section 3.4). Disc spectrum calculated 
under the assumption that the local emission is well approx- 
imated by black body is shown in Fig. 8. More advanced 
disc spectra are also available in the literature (e.g. Czerny 
& Elvis 1987, Laor & Netzer 1989, Ross, Fabian & Mineshige 
1992, Shimura & Takahara 1995) but as the accuracy of the 
description is still problematic due to simplifications in the 
description of the disc vertical structure and the radiative 
transfer, we decided to use the simplest model for the pur- 
pose of illustration. 

The examples of the resulting spectra (without Comp- 
tonization by the corona) are shown in Fig. 8. They were 
calculating assuming Kcplcrian disc locally radiating as a 
black body but both the mass loss (see Section 3.4) and the 
decrease of the internal dissipation due to the presence of 
the corona (equation 18) were taken into account. 

As the accretion rates drops by a significant factor on 
the way to the black hole the overall spectrum is flatter in 
the optical band than the usual spectrum of a stationary 
conservative Keplerian disc. 

In our case the mass loss (and subsequently the energy 
index a_B, Fe ~ E^"^ in UV) is larger for larger accretion 
rate. However, this trend may reverse if there is any change 
of the viscosity parameter with the corona temperature or 
bolometric luminosity. 

3.6.3 X-ray spectrum 

For most values of external accretion rates we predict sig- 
nificant Comptonization in the corona if the mass loss from 
the disc is taken into account. The values of the Comp- 
ton parameter y are typically 0.5 — 0.7, as can be seen from 
Fig. 6 (2/ ~ 1 — ^). As our computations have still approx- 
imate character we did not compute the resulting spectra 
taking into account the radial properties of the corona. In- 
stead, we calculated the accretion rate at 10 ff.schw appro- 
priate for assumed external accretion rate, computed the 
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Figure 8. Two examples of the spectral shape of the opti- 
caUy thick, local black body emission of a stationary, Keple- 
rian accretion disc with both the effect of g < 1 and decreas- 
ing (with r decreasing) local accretion rate due to the outflow 
(solid curves). The external accretion rates are 0.001 and 0.01 
in units of lO^^g/s. The dashed curves show effects of Comp- 
tonization of the spectra by the corresponding coronae (see text 
for details). For comparison, the spectrum of a conservative disc 
(M(r) = const) without a corona = 1) is shown with a dotted 
curve for the same two values of the accretion rate. 

corona parameter at this radius and Comptonized the total 
disc spectrum assuming a uniform corona with properties as 
discussed above. At the present stage we do not include the 
X-ray radiation created in the corona but reflected by the 
disc before reaching an observer. 

The result of the Comptonization is the formation of 
the hard X-ray power law tail with energy index 0.92 for 
lower value of the external accretion rate (0.001 x lO^^g/s) 
and 1.35 for the 10 times higher value. This first value 
corresponds well to the observed spectra of Seyfert galaxies 
if they corrected for the reflection component (Pounds et 
al. 1990; Nandra &; Pounds 1994). The second value may 
be adequate for a few exceptional Seyfert galaxies like RE 
1034-1-39 (Pounds et al. 1995) and some quasars. 

The predicted spectra, however, depend strongly on the 
estimates of the mass loss from the disc surface which in 
turn strongly depends on the viscosity parameter (see Sec- 
tion 3.4). If the actual mass loss is lower the corona is weaker 
and the hard X-ray power law of the required slope would 
not form by Comptonization in a continuous corona. In 
that case the hard X-ray emission should originate in com- 
pact active regions developing in the corona most probably 
due to reconnections of the magnetic field lines, as recently 
discussed by Haardt, Maraschi & Gliisellini 1994 and Stern 
et al. 1995. 

In any case the Comptonization by the corona leads 
to modification of the high frequency tail of the big bump 
leading to a power law shape instead of thermal exponen- 
tial cut-off. Actually, in most sources with very strong soft 
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X-ray excesses the spectrum looks like an extension of the 
thermal emission of the big bump but the slope is a power 
low instead of an exponential cut-off (e.g. Czerny & Elvis 
1987, Walter & Fink 1994). On the basis of EUVE and 
ROSAT WFC observations Marshall, FVuscione & Carone 
(1995) argue that the typical Seyfert spectrum in the X and 
UV range is given approximately by a power law with en- 
ergy index ~ 1.16 — 2.07. Such a shape of the spectrum is 
easily explained by Comptonization by the plasma with a 
low to moderate value of the Compton parameter. Our con- 
tinuous corona is therefore in a natural way identified with 
this medium. 

EYirther more detailed modeling is necessary before we 
can distinguish between different scenarios of the hard X- 
ray production and the origin of soft X-ray emission for a 
given source. Variability studies will be particularly helpful 
for that purpose (e.g. Czerny, Zycki &; Loska 1995). 



4 RESULTS FOR THE CORONA WITH 
RADIAL ADVECTION 

4.1 Role of advection 

The set of our equations describing the heating/cooling of 
the corona (see Section 2.3) does not take into account the 
radial energy transport. The most important term is the 
radial advection; its importance in the energy budget of the 
fast accretion flow was discussed in a number of papers (Mu- 
chotrzeb k, Paczynski 1982; Muchotrzeb 1983, Taam & Lin 
1984; Abramowicz ct al. 1988; Chen & Taam 1993; Narayan 
& Popham 1993; Narayan & Yi 1994,; Chen 1995). In the 
case of coronal accretion the corresponding expression for 
advective heating is given by 



47rr_Fa. 



(1 



(28) 



where (1 — QM is the local accretion rate in the corona, Ti 
is the temperature (i.e. the ion temperature in the case of 
two temperature medium and S is the mean entropy in the 
vertical direction. For a gas pressure dominated medium 
this equation can be conveniently expressed as 

.dlnTi dlnP, 



47rr^i^adv = (1 - OMcib- 



(2.5- 



(29) 



rf In r d In r 

The sign of this term depends on the enthropy gradient in 
the flow. We can estimate it assuming that the advection is 
negligible, i.e. on the basis of solutions presented in Section 
3. The value of the coefficient 6 is approximately 



-3/4, 



(30) 



mostly due to the fact that the ion temperature is a de- 
creasing function of radius and the cooler gas flows into the 
hotter coronal region thus acting as a cooling mechanism. 
We adopt the expression (29) with the value of S given by 
(30) and we recalculate the structure of the corona adding 
the advection term to the viscous dissipation term in equa- 
tion (16). Similar method was used while computing the 
role of advection in the disc itself (Chen 1995). 

4.2 Properties of the corona 

The inclusion of the advection term, as described in 
Section 4.1., does not change the results significantly. The 
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Figure 9. We illustrate the importance of the radial advection 
in our model showing the fraction of the energy liberated in the 
disc as a function of accretion rate without advection (continuous 
lines) and with advection term included (dashed lines); the corona 
dissipates a fraction 1 — 5- Curves are labelled by the viscosity 
parameter a. For better clarity, only two cases with advection 
are shown (a — 0.1 and 1.0). Other parameters: M — 10* Mq, 
r = lOiischw 

advection is unimportant at low accretion rates. Increas- 
ing accretion rate we still first meet the topological prob- 
lems with solutions - first the transonic curve does not pass 
through the first singular point and soon after the two sonic 
points converge (see Section 3.f) - and the solutions cease 
to exist before the advection term has a chance to dominate 
the thermal balance of the flow. 

Therefore the only result of the advection term is to 
extend the solutions slightly towards higher accretion rates 
and to lower the fraction of the energy liberated in the disc 
for a given accretion rate, but the differences are practically 
negligible. To allow the comparison with the previous case 
wc show the efficiency plot for two values of the viscosity 
parameter a (Fig. 9). Other properties of the corona also 
do not practically differ from those computed neglecting the 
advection term (see Fig. 4) so we do not repeat the plots in 
this Section. 

Our results strongly differ from advection-dominated 
solutions for the entire disc. This is due to the fact that 
stratified flow is cooled far more efficiently than an unstrati- 
fied one as the optically thick layer supplies copious photons 
to the corona. 



5 DISCUSSION AND CONCLUSIONS 
5.1 Luminosity states 

The dependence of the disc/corona solution on accretion 
rate at a given radius show three distinct states. At very low 
accretion rates the corona does not develop, at intermediate 

accretion rates a strong corona forms and at high accretion 
rates the corona becomes relatively weak. Such states may 
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correspond to observed states: very low state, low state and 
high state. 

Very low state was observed in the case of a galactic 
source V404 Cyg and suggested in the case of accretion 
at the central black hole of the Galaxy (Falcke & Heinrich 
1994). Such a state is characterized by very soft thermal 
emission, possibly corresponding to standard accretion disc 
although the typical temperature is somewhat higher than 
expected on the basis of accretion rate appropriate for the 
observed luminosity. 

Low state is well observed in a number of black hole 
candidate galactic sources and show general similarity with 
the spectra of most Seyfert galaxies (see Mushotzky et al. 
1993). Its characteristic property is the presence of hard 
X-ray power law with the energy index ~ 0.9 (or sometimes 
flatter) extending to a few hundreds of keV. This compo- 
nent contains considerable or even dominating fraction of 
the total bolometric luminosity of the source. 

High states are also well defined in the case of black 
hole galactic candidates. The spectrum is again relatively 
soft, most probably due to the thermal emission of optically 
thick medium and it is approximately modeled by accre- 
tion disc although the temperatures are somewhat higher 
than predicted by standard accretion discs due to the Comp- 
tonization inside optically thick discs. Similar spectra are 
characteristic for quasars and exceptional Seyfert galaxies 
(e.g. RE 1034-1-39, Puchnarewicz et al. 1995; Pounds, Done 
and Osborne 1995). 

The coincidence of these trends with the predictions 
of our model is encouraging although several simplifications 
were made in our computations so it is difficult to expect de- 
tail quantitative agreement between the data and the model 
a this stage. 

5.2 Spectra 

The optical/UV/X-ray spectra computed on the basis of our 
model are in agreement with the trends observed in the data. 
The optical spectra are softer than predicted by standard 
accretion disc model due to the significant mass loss from the 
disc surface and subsequent decrease of the accretion rate 
towards the central black hole. The slope of the hard X-ray 
power low is also reproduced for external values of accretion 
rate about 0.01 of the Eddington value due to formation of 
the strong corona. 

However, our results should not be treated as a proof 
that the accreting continuous corona is actually the best ex- 
planation of the origin of the X-ray emission. Although the 
model predictions at a given radius do not crucially depend 
on the assumptions, including the value of the viscosity pa- 
rameter, the predictions of the radial dependencies, partic- 
ularly of the mass loss, do depend strongly on the viscosity. 
Further work is needed before any distinction between ac- 
creting continuous and patchy magnetically driven corona 
(e.g. Haardt et al. 1994, Stern et al. 1995) could be made. 

5.3 Outflow 

In our model the accretion flow into the black hole is ac- 
companied by very strong outflow. Such a significant out- 
flow of the material seems to be one of the observed basic 
characteristics of nuclear activity. Radio-loud objects have 



collimated outflows whilst similar but uncoUimated outflow 
is argued for in radio-quiet objects (e.g. Stocke et al. 1992), 
both quasars and Seyfert galaxies. 

There are several direct evidences of the outflow from 
the central regions of active galaxies although the local- 
ization of the origin of the outflow is hardly constrained. 
Absorption features due to the partially ionized medium 
present along the line of sight to the nucleus (i.e. warm 
absorber) detected in a number of Seyfert 1 and Narrow 
Line Seyfert galaxies (e.g. Halpern 1984; Nandra & Pounds 
1994; Fiore et al. 1993; Turner et al. 1993) are most proba- 
bly due to the outflowing gas (Reynolds & Fabian 1995). 
Similar features are observed in BL Lac objects (Made- 
jski et al. 1991). ASCA observations of NGC 4051 (Mi- 
hara et al. 1994) suggest the velocity ~ 10000 km/s and 
lower values (~ 1000 km/s) were determined for NGC 4151 
(Kriss et al. 1992). Evidences of an outflow with velocities 
~ 3500 - 5000 km/s are also seen in Seyferts PG 1351-1-64, 
Ton 951 and III Zw 2 (Stocke et al. 1994). The significant 
change of absorption features in three weeks observed for 
MCG-6-30-15 (Fabian et al. 1994) indicate that the warm 
absorber is located within a few light weeks. It may have the 
covering factor close to unity (e.g. George, Turner &; Netzer 
1995 for NGC 3783). Further out the outflow may continue 
through the BLR (although in the case of the wind from the 
inner edge of the dusty /molecular torus the outflow is not a 
continuation of the nuclear wind, Krolik &; Begelman 1986, 
Balsara & Krolik 1995) to NLR (e.g. Krolik & Vrtilek 1984; 
see also Wilson 1994, Mulchaey, Tsvetanov & Wilson 1994). 
The mass flux in that region was estimated as ~ IMo/yr 
(Krolik & Vrtilek 1984). 

Indirect support for the nuclear wind comes from the 
Broad Line Region model developed by Smith & Raine 
(1985) and Mobaslier & Raine (1989). Clouds in their model 
forms due to the collision of the nuclear wind with another 
wind developed from outer parts of the disc. They expect 
nuclear wind velocities of order of 10000 — 50000 km/s on the 
basis of line profiles and intensity ratios. The velocities of 
the outflow of the nuclear wind predicted by our model are 
also in that range which is encouraging for future studies. 
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APPENDIX A 

We show here how equations (9) and (10) can be partly 
solved analytically. Surprisingly the second differential 
equation (10) can be integrated without knowing the so- 
lution for rh explicitly. We only need to multiply the second 
differential equation by an integrating factor which is 
in this case. 

We now obtain an implicit function in P for the solution 
of the differential equation: 

F(.,P) = C/i+ "^'^-^ 

(Al) 





RschwC^rnii 



2kT 



Vr''' + z 



ckT ^ 
' + in P : 



where C\ is the integrating constant. If the solution shall 
go through the sonic point we obtain for this case 

m^fcT -Fsoit'>'es?nii , 



F{z,P) 



2 2P2,„jj 



RschvC mil 



ckT 



-{z 



tlic) 



2kT 



In 



V J^sonic / 



■- 0. 



\/r'^ -\- z. 



2 

sonic 



+ (A2) 



To obtain a solution for m and subsequently for P we 
assume first that the viscosity a = 0, i.e. Vt = 0. In this case 
the two differential equations decouple and we can solve for 
m. For the solution we obtain 

C2 



r2 + z^ 



(A3) 



where C2 is the second integrating constant. To solve the 
differential equations for a ^ we take the ansatz 



= 2 I 2 ' 
y.2 _|_ ^2 

which yields 

df _ 2vrrmn 
dz 



kT 



iz) 



kT / ^^^''^^ 



(A4) 



(A5) 



The first term in this equation represents the mass flux from 
the disc to the corona whilst the second term gives the net 
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effect of the radial flow in the corona itself. 

In principle we can insert the complete solution for m 
into equation (Al) which yields us an equation only for P. 
This procedure yields us an integral equation for P which 
can be transformed into an implicit differential equation only 
in P. 

Since the solutions for P{z) for different viscosity pa- 
rameters a do not differ very much from the solution for 
Q = we can interpret a as a small distortion and try to 
obtain an approximate solution by Taylor expansion. 

For low accretion rates (M < 0.1 x lO^^g/s) we can 
neglect the term rh^kT/2miiP^ in equation (A2) for small z 
and obtain a very good approximation for the pressure term 

Fsott Ke 



P{z) =Psomc exp + 



cT 



c) 



+ 



flschwC^"^H 

2kf 



Vr^ + z^ y^r2 + zl 



(A6) 



for zo < z ^ Zsonic- The error at the location zo is usually 
much smaller than 1%. The error increases with larger z and 
at the location «sonic in equation (A6) the pressure is over- 
estimated by a factor e^^'^. However, for very low accretion 
rates (M < 0.01 x lO^*' g/s) equation (A6) yields excellent 
results because the relevant range for the pressure terms is 
in the range zo < z ^ «sonic, where P is sufficient large. 

In general we can state that the pressure decreases ex- 
ponentially with larger z in the range zo ^ z ^ ^sonic. Con- 
sidering equation (A4) and equation (A5), the same holds 
for m for a > 0. 



APPENDIX B 

In Appendix A we learned that we can decouple the two 
differential equations (9) and (10) for rn{z) and P{z) and 
obtain a differential equation only in m or P. 

To compute a contour plot for P{z) it is not at all clear 
that we can choose the second integrating constant C2 in 
equation (A5) independently of the integrating constant Ci 
in equation (Al). 

We derive here the curve where dP/dz — » 00 in the 
contour plots in Figure 1. This curve fully determines the 
initial values for P and rh of the contour plot if we determine 
the transonic solution. To obtain dPjdz — > cxs in the contour 
plot we have the condition P = m^JkT /jnu. Inserting now 
the solution of m of equations (A4) and (A5) yields 



(r^ + z^)Pl 




Ci -I- 



P{z)dz 



kT/rriH 

^0 

After differentiating equation (Bl) we obtain 

-.P. 



dP , 2 2n „ Sii^r 
— {r^+z^) + 2zP: 



dz 



(Bl) 



(B2) 



mii 



For this case the variables can be easily separated and the 
differential equation can be solved. Finally we obtain the 

solution, 



J'initial(a) 



+ z- 



exp 



2Vr 



z ' 

arctan(-) 
/ r 



(B3) 



The integrating constant C3 has to be calibrated in such 
a way that Pinitiai('Zsonic) = fsonic. The initial values for 



m ca n be de termined subsequently by the condition rh 
P/y/kT/mn. 



APPENDIX C 

We can derive simple analytical expressions for the param- 
eters of the corona as a function of a mass of the central 
black hole, M, disc radius x measured in Scliwarzscliild units 
{x = r/Rschw) and accretion rate in units of the Eddington 
accretion rate with efficiency 1/16. For that purpose we as- 
sume that the dynamical terms and the advection term are 
negligible, the ion temperature is much larger than the elec- 
tron temperature, the corona structure ran be approximated 
using vertically averaged quantities and the nonlinear terms 
in energy balance equations are negligible. We also approx- 
imate ln((.3 — 0/(1 + 0) with 1 — ^ which is accurate for 
^ ~ 1). The boundary condition is included through /(r) 
(see equation 14a). 

Under these assumptions the parameters of the corona 

are: 



1 - ^ = 0.017a-i/VW"'^'^"'^ 

Ti = 5.2 X 10''a-''^f{r)"^m''\-'''', 

Te = 1.3 X 10'a-''^f{r)-''^m-''\''\ 

and 

T = 0.18. 



(CI) 
(C2) 
(C3) 
(C4) 

(C5) 



As we see, the properties of the corona do not depend on 
the mass of the central object so the same model can be 
applied to AGN and galactic black hole candidates. This is 
in contrast with the cold disc properties like temperature or 
the gas pressure to the total pressure ratio. 



APPENDIX D 

In this Appendix we derive a set of equations which in first 

order represents the global solutions for the Figs. 3 and 4. 
The equations are needed to obtain the results for different 
r with higher accuracy that provided by formulae given in 
Appendix (C). This higher accuracy is essential to calculate 
the vertical vs. radial flow (see equation 26) for different 
radii. Since equation (A6) represents a very good solution 
for the pressure for small z we have P{zo) ~ P{0) = Po and 
we can write for z <C Zsonic 



which yields 



^{z) = Po exp 



MGrnu a 
kTi 2^ 



(Dl) 



In addition we neglected the term with Fgoft in equation 
(Dl). Equation (Dl) is a good approximation for small z, 
but a bad one for large z. However, the overall integral 
over the pressure should represent a good approximation 
because the main quantity of the integral is contributed by 
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small z. Therefore the decrease of the function for large z is 
not important for this purpose. Using now the relation 



oo 

/ 



i_{z)dz 




(D2) 



and the relation 



oo 

/ 



'approx(^)'^^ — 2 




yields us for equation (16) 

(1 - O-Pto. 2 — "V 2 
and for equations (20) and (21) 

3 - ^ _ f¥ Po Kes 4fcTe 
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1 + ^ V 2 ^kTilmn "^<=c' 
Further we obtain for equation (17) 



(1 - i)Fu 
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X 2.44 X 10^' In A 
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VmeC^ 
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(D3) 



(D4) 



(D5) 



(D6) 



Finally we need a fourth equation, which is represented by 
equation (22) and which can be taken without alteration. 



Ho 



ViosKy 



2cPo 



(D7) 



Equations (D4)-(D7) represents the set of equations which 
yields us the global solution for the set of parameters M, M, 
a, r. The four equations can be solved easily numerically 
and they yield us the four unknown quantities ^, Ti, Te and 
Po- Also equations (D4)-(D7) can be used to derive the 
equations in Appendix C. For example, equations (D4) and 
(D7) yield immediately that Ti ex in first order. 
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